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Prufungsantrag gem. § 44 PatG ist gestellt 

@ Verfahren zur Berechnung von Strukturen mit Hilfe der Finiten-Elemente-Methode 

@ Verfahren zur Berechnung von Strukturen mit Hilfe der 
Finiten-Elemente-Methode mit den Schritten: Idealisie- 
rung der zu berechnenden Struktur durch Volumenele- 
mente und Festlegung der Eck-Knotenpunkte in der idea- 
lisierten Struktur, Zerlegung der Volumenelemente mit 
Hexaeder-, Pentaeder- und Pyramiden-Form in Tetraeder, 
so dafc die zu idealisierende Struktur eine Teilelement- 
Struktur aus Tetraedern ist, Beschreibung des Verschie- 
bungszustands f u r jedes Tetraeder-Element durch ein 
vollstandiges Polynom dritter Ordnung je Koordinaten- 
richtung, Bestimmung der Polynom-Koeffizienten durch 
Zuordnungen zwischen den Knotenfreiheitsgraden, den 
Faktoren zur Beschreibung der individuellen Tetraeder- 
■ geometrie und den jeweiligen Polynom-Koeffizienten, 
wobei als Knotenfreiheitsgrad fur jeden Knoten-Eckpunkt 
der Verschiebungsvektor und die drei Ableitungen dieses 
Vektors und der Verschiebungsvektor fur die Flachenmit- 
ten-Knoten jedes Tetraederelements verwendet wird, Eli- 
mination der Flachenmitten-Knotenpunkte, Integration 
und Losung des gestellten mechanischen Problems mit- 
tels der Systemmatrix. 
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Beschreibung 

Die Erfindung betrifft ein Verfahren zur Berechnung von Slrukturen mil Hilfe der Finiten-Elemente-Methode, das sich 
insbesondere auch zur Berechnung von Schalenstrukturen eignet. 

5 Finite Elemcnteverfahren laufcn ublichcrwcisc in folgenden drei Schrittcn ab: cinem Aufbereitungsschritt zur Idcaii- 
sierung der zu berechnenden Struktur, einem Analyseschritt, d. h. einem Berechnungsschritt zur Berechnung der in der 
zu berechnenden Struktur auftretenden Verformungen, Spannungen, Dehnungen, StabilitatsgrbBen und dynamise hen La- 
sten, und einem Auswertungsschritt zur Darstellung der Ergebnisse. In dem Berechnungsschritt werden bei der Finite- 
Elemente-Methode allgemein auf der Basis einer Vielzahl von finiten Elementen, die im Vorverarbeitungsschritt zur 

10 Tdealisierung der zu berechnenden Struktur in einer raumlichen Anordnung definiert worden sind, Systemmatrizen und 
zusammen mit den auBeren Kraften Lastmatrizen gebildet, urn daraus VerschiebunggroBen an den Knotenpunkten der 
gewahlten finiten Elernente zu berechnen. Die Systemmatrizen werden dabei aus sogenannten Elementmatrizen zusam- 
mengesetzt, d. h. aus Matrizen, die aus jedem einzelnen zur Idealisierung der Struktur verwendeten finiten Element er- 
mittelt werden. 

15 Zur Bildung der Elementmatrizen wurden bisher Volumenclemcnte rnit Polynom-Ansatzcn niedriger Ordnung, d. h, 
bis zur zweiten Ordnung, verwendet. Dadurch konnte die berechnungstechnische Funktionsfahigkeit bei der Losung der 
Ansatze beherrscht werden. Bei diesem Verfahren hatte der Benutzer die kinematischen Freiheitsgrade fur jeden System- 
knoten der idealisierten Struktur manuel, dabei individuell festzulegen und zu steuern. 

Bei dem Verfahren nach dem Stand der Technik werden ftir die finiten Elernente unterschiedliche geometrische For- 

20 men und kinematische Formulierungen verwendet. Als geometrische Formen konnen dabei bei spiels weise Schcibenele- 
mente, Plattenelemente, Schalenelemente, Stabc oder Balken oder auch Volumenelemente vorgesehen sein. Die Aus- 
wahl der entsprechenden geometrischen Formen fiir die finiten Elernente erfolgt in Abhangigkeit der Gestalt und mecha- 
nischen Funktion der zu idealisierenden Struktur bzw. deren Strukturteile. Dabei werden beispielsweise bei schlanken 
Strukturteilen als finite Elernente Balken oder Stabe verwendet, wahrend bei flachigen Strukturen bzw. Strukturteilen 

25 Scheiben, Platten oder Schalen verwendet werden. Die Auswahl der jeweiligen geometrischen Formen steht im Ermes- 
sen des Benutzers, d. h. sie bemiBt sich nach der Erfahrung und Einschatzung des Benutzers. In vielen Fallen stellt sich 
erst nach Durchfuhrung eines Berechnungsdurchlaufes und nach Analyse der Ergebnisse sowie nach Vergleich mit Ver- 
suchsergebnissen heraus, welche Qualitatdie Idealisierung der Struktur und deren Berechnung aurweist. In vielen Fallen 
ist eine nochmalige verbesserte Idealisierung der Struktur nbtig, die jedoch wiederum manuell aufgrund qualitativer Be- 

30 wertungen der Berechnungs- und Versuchsergebnisse zu erfolgen hat. Diese Verfahren nach dem Stand der Technik sind 
somit sehr arbeitsaufwendig und machen umfangreiche Versuche an der zu berechnenden Struktur erforderlich, urn die 
Berechnungsergebnisse zu perfektionieren. Zudem bieibt die Qualitat der insgesamt erreichbaren Ergebnisse begrenzt 
und benutzerabhangig. 

Vor Anwendung von Finite-Elemente-Methoden wurde zur Berechnung von Strukturen eine sogenannte KraftgroBen- 

35 methode angewendet, bei der als Unbekannte aus auBeren Kraften abgeleitete innere Krafte verwendet werden. Aus die- 
sen inneren Kraften wurden die Spannungen, Dehnungen und Verschiebungen ermittelt, wobei die Verschiebungen nur 
ungenau ermittelt werden konnten, da diese nur indirekt errechenbare GroBen sind. Mit einer Verbesserung der Lei- 
stungsfahigkeit der verwendeten Rechenanlagen wurde die KraftgrbBenmethode durch die Finite-Elemente- Methode er- 
setzt, wobei eine sogenannte Verschiebungsmethode angewendet wurde, bei der Verschiebungen als Unbekannte ver- 

40 wendet werden. Bei dieser Methode werden in Abhangigkeit der ermittelten Verschiebungen die gesuchten Dehnungen 
und Spannungen berechnet. AuBerdem sind nach dem Stand der Technik auch Kombinationen aus der KraftgroBenme- 
thode und der Verschiebungsmethode bekannt. Derartige kombinierte Verfahren bieten jedoch keine signifikante Verbes- 
serung hinsichtlich der Ergebnisgenauigkeit im Vergleich zur reinen Verschiebungsmethode. 
Ein Nachteil sowohl der Verschiebungsmethode als auch der kombinierten Verfahren ist, daB sie einen groBen Rechen- 

45 aufwand benotigen. Ein weiterer Nachteil ist, daB sie zu numerischer Empfindlichkeit neigen, d. h. die Ergebnisqualitat 
hangt im hohen MaBe von der richtig gewahlten Feinheit der Idealisierung der zu berechnenden Struktur ab, wobei sich 
diese nach der Einschatzung des Benutzers bestimmt. Bei sehr feinmaschig idealisierten Strukturen werden haufig auch 
unverhaltnismaBig hohe Verformungswerte errechnet. Mit Verfahren nach dem Stand der Technik konnten weiterhin 
viele Stabiiitatsprobleme z. B, Beulprobleme nicht gelbst werden. Insbesondere die Berechnung von gekrummten und 

50 gefalteten Schalen konnte in vielen Fallen nicht gelbst werden, da die Faltung der Struktur zu Kompatibilitatsdefekten 
bei der Analyse fiihrte, 

Es ist daher die Aufgabe der Erfindung, ein Verfahren zur Berechnung von Strukturen mit Hilfe der Finite- Elemente- 
Methode zu schaffen, die bei im Vergleich zum Stand der Technik vergleichbarer Rechenkapazitat zu einer Steigerung 
der Ergebnisqualitat fuhrt, wobei die numerische Stabilitat insbesondere bei schlanken Elementen gewahrleistet sein 
55 muB. Mit dem zu schaffenden Verfahren sollen auch mathematisch aufwendige mechanise he Stabiiitatsprobleme gelbst 
werden konnen. AuBerdem soli der EinfluB der benutzerabhangigen \brgehensweise erheblich reduziert und eingegrenzt 
werden. 

Diese Aufgabe wird mit den Merkmalen des Anspruchs 1 gelost. Alternative Ausfiihrungsformen sind in den Unter- 
anspruchen angegeben. 

60 Nach dem erfindungsgemaBen Verfahren entfallt die Auswahl der Element-'iypen, d. h. der geometrischen Formen der 
finiten Elernente, bei der Idealisierung der zu berechnenden Struktur, da erfindungsgemaB ausschlieBlich Volumenele- 
mente zur Idealisierung der Struktur verwendet werden. Das erfindungsgemaBe Verfahren betrachtet die zur Idealisie- 
rung definierten Volumenelemente, soweit sie nicht als Tetraeder definiert worden sind, als Verbun del erne nte und setzt 
sie aus Tetraedern zusammen. Bei schlanken Strukturen und Schalenstrukturen ist es ausreichend, uber die Dicke der 

65 Struktur nur eine Lage von Volumenelementen vorzusehen. ErfindungsgemaB werden also keine ein- und zweidimensio- 
nalen Elernente, also auch keine Schalenelemente oder Stabelemente, zur Idealisierung der Struktur verwendet. Dadurch 
wird die Durchfuhrung der Idealisierung der zu berechnenden Struktur systematisiert und vereinfacht, ohne daB die Er- 
gebnisgenauigkeit dadurch beeintrachtigt wird. Durch die Verwendung ausschlieBlich von Volumenelementen als Finite 
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Elemente fur die idealisierte Struktur werden die nach dem Stand der Technik ub lie hen Einschrankungen der Kirch- 
hoff schen Biegung sowie der sogenannten Schubbiegung nach Reissner-Mindlin aufgegeben. Nach der Kirch hoffschen 
Vorgehensweise wird angenommen, daB der ebene, normale Querschnitt der idealisierten Struktur bei einer Beanspru- 
chungen, z. B. einer Biegung, normal und eben erhalten bleibt. Nach der Reissner-Mindlin- Annahme werden Schubver- 
zerrungcn bei Biegcbeanspruchung zwar beriicksichtigt, aber normale, ebene Qucrschnittc bleiben auch erhalten. Durch 5 
Verzicht auf diese Annahmen oder Einschrankungen zusammen mit der Verwendung von Polynom-Ansatzfunktionen 
dritter Ordnung fiir den Verschiebungszustand jedes Tetraeder-Teilelements bezUglich jeder Koordinatenrichtung lassen 
sich Unstetigkeiten in der Struktur berechnungstechnisch genauer modellieren. Nach dem Stand der Technik ist die An- 
nahme, daB der Querschnitt der zu berechnenden Struktur eben bleibt, eine Einschrankung, durch die die Abweichungen 
bei der Berechnung der Strukturen groB werden, was zu unkalkulierbaren Ungenauigkeiten bei der Berechnung der 10 
Struktur fiihren kann. Auch bei der Beschreibung des Verschiebungszustandes jedes Tetraeder- Elements, wie es in einem 
nachfolgendem Schritt erfolgt, werden diese Einschrankungen auch nicht - wie im Stand der Technik - indirekt wieder 
eingefuhrt, da hierfiir erfindungsgemaB Polynom-Ansatzfunktionen dritter Ordnung verwendet werden, die geeignet 
sind, den allgemeinen Verschiebungszustand ohne Einschrankungen zu beschreiben. 

Wciterhin werden erfindungsgemaB fiir die Verformungcn Polynomansatze verwendet, die zumindest dritte Ordnung 15 
aufweisen. Dadurch ergibt sich zum Vergleich beim Stand der Technik eine wesentlich bessere Annaherung der realen 
Verformungszustande fur die zu berechnende Struktur. 

ErfindungsgemaB werden fur die Eck-Knotenpunkte der idealisierten Struktur als Freiheitsgrade die TYanslationen in 
alien drei Richtungen des Raumes sowie die Ableitungen erster Ordnung zu diesen Trans lationen vorgesehen. Zusatzlich 
werden nach der Erfindung die Verschiebungen der Flachen-Mittenknoten der Tetraeder-Idcalisierungs- und -Tcilemcnte 20 
vorubergehend eingefuhrt. Demgegeniiber sind die Freiheitsgrade bei Verfahren nach dem Stand der Technik abhangig 
von den jeweils verwendeten Finite-Elementtypen zu definieren. Aus diesem Grund muBten bei Verfahren nach dem 
Stand der Technik die Freiheitsgrade durch den Benutzer mit relativ groBen tlberlegungsaufwand eingegeben oder ge- 
steuert werden, da beispielsweise an Ubergangsstellen verschiedenartiger Finite-Element-iypen, z. B. beim Ubergang 
von Stabelementen zu Scheibenelementen, schwer zu ermittelnde Ubergangs-Randbedingungen vorzugeben sind. Ann- 25 
liche Problemc ergeben sich bei Verwendung von Elementtypen, die in bestimmten Richtungen nur eine sehr geringe 
oder gar keine Steifigkeit aufweisen. Bei dem Verfahren nach der Erfindung dagegen sind wegen der Verwendung von 
Volumenelementen alle Steifigkeiten in alien Richtungen a priori vorhanden, so daB keine von den verwendeten Ele- 
menttypen abhangige Uberlegungen anzustellen sind. Vielmehr ergibt sich die Moglichkeit einer einheitlichen Vorge- 
hensweise fiir samtliche Strukturtypen, die auch automatisiert werden kann. 30 

Ein weiterer Vorteil der Erfindung ist, daB durch die Verwendung ausschlieBlich von Volumenelementen mit heute ub- 
lichen Konstruktions-Werkzeugen auf die Ergebnisse des erfindungsgemaBen Verfahrens angewendet werden konnen, so 
daB das im erfindungsgemaBen Verfahren entwickelte dreidimensionale Modell natiirlich, d. h. formgetreu ubemommen 
werden kann. 

Im folgenden wird die Erfindung anhand der beigeftigten Figuren beschrieben. Es zeigen: 35 

Fig. 1 eine schematische Darstellung eines Hexaeders als eines erfindungsgemaB zur Struktur-Tdealisierung verwen- 
deten Verbund- oder Volumenelements in perspektivischer Darstellung, das nach einer ersten Variante in funf Tetraeder - 
vier peripheren und einem zentralen Tetraeder im Volumeninnern - nach dem erfindungsgemaBen Verfahren geteilt wird, 
wobei der fiinfte, vollstandig im Inneren des Hexaeders liegende Tetraeder separat dargestellt ist, 

Fig. 2 eine perspektivische Darstellung eines Hexaeders als eines erfindungsgemaB zur Struktur-Idealisierung verwen- 40 
deten Volumenelements, das nach einer zweiten Variante - in systematischer Erganzung/Abwandlung der ersten Variante 
nach der Fig. 1 - nach der Erfindung in funf Tetraeder geteilt wird, wobei der fiinfte, vollstandig im Inneren des Hexa- 
eders liegende Tetraeder separat dargestellt ist, 

Fig. 3 eine perspektivische Darstellung eines Pentaeders, der alternativ als Verbundelement vorgesehen werden kann, 
der erfindungsgemaB nach einer ersten Variante in drei Tetraeder geteilt wird, 45 

Fig. 4 eine perspektivische Darstellung eines als Verbunedelement verwendeten Pentaeders, das neben der Teilungs- 
naht der ersten Variante der Fig. 3 nach einer zweiten Variante aus drei entegegengesetzt orientierten Tretraedem geteilt 
wird, 

Fig. 5 eine perspektivische Darstellung eines als Verbundelement verwendeten Pyramide, die erfindungsgemaB nach 
einer ersten Variante in zwei Tetraeder geteilt wird, 50 

Fig. 6 eine perspektivische Darstellung einer als Verbundelement verwendbaren Pyramide, die erfindungsgemaB nach 
einer zweiten Variante aus zwei systematise h altemativen Tetraedern zusammengesetzt ist, 

Fig. 7 eine Draufsicht auf einen Tretraeder, der erfindungsgemaB sowohl als Verbundelement als auch als Teilelement 
eines Verbundeiements nach den Fig. 1 bis 6 verwendbar ist und 

Fig. 8 eine schematische Darstellung des Ablaufs des erfindungsgemaBen Verfahrens. 55 

Die Fig. 1 bis 7 zeigen verschiedene Ausfuhrungsformen fiir Volumenelemente, die dem erfindungsgemaBen Verfah- 
ren zur Berechnung von Strukturen mittels der Finite-Element- Methode verwendbar sind. Dabei werden in der zu idea- 
lisierenden Struktur Knotenpunkte definiert, die die Eckpunkte fur die jeweils verwendeten \fclumenelemcnte bilden. In 
dem erfindungsgemaBen Verfahren werden die als Idealisierungs-Elemente verwendeten Volumenelemente, soweit diese 
nicht selbst Tetraeder sind, in Gruppen aus Tetraedern zerlegt. Da diese Volumenelemente aus Tetraedern zusammenge- 60 
setzt sind, werden sie als Verbundelemte bezeichnet. In den Fig. 1 und 2 ist beispielhaft dargestellt, wie ein als Verbund- 
oder Volumenelement verwendeter Hexaeder aus fiinf Tetraedern aufgebaut sein kann. Der Aufbau aus Tetraedern ist in 
den Fig. 3 und 4 fiir den Pentaeder und in den Fig. 5 und 6 fur die Pyramide gezeigt. Das erfindungsgemaBe finite Ele- 
ment- Verfahren erfolgt dann auf der Basis der Tetraeder- Anordnung en. In der Darstellung der Fig. 1 bis 7 stellen die dik- 
ken durchgezogenen Linien sichtbare Kanten, die dunnen durchgezogenen Linien sichtbare Diagonalen, die gestrichel- 65 
ten Linien unsichtbare Kanten und die gepunkteten Linien unsichtbare Diagonalen dar. 

Das in der Fig. 1 dargestellte Idealisierungs-Element ist ein Hexaeder 1. In der Fig. 1 ist dargestellt, wie dieses erfin- 
dungsgemaB in Tetraeder zerlegt wird. Dabei bleiben die Eckpunkte 2, 3, 4, 5, 6, 7, 8, 9 als Eck-Knotenpunkte fur die Te- 
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traeder-Teilelemente 11, 12, 13, 14, 15 erhalten. 

In der Fig, 2 ist eine alternative Zerlegung eines Hexaeders 1 als Idealisierungs-Element dargestellt. Auch hier bleiben 
die Eck-Knotenpunkte 2, 3, 4, 5, 6, 7, 8, 9 bei der Zerlegung des Hexaeders 1 in die Tetraeder 11, 12, 13, 14 erhalten. 

Eine Zerlegung eines Pentaeders 21 mit den Eckpunkten 22, 23, 24, 25, 26, 27 in die Tetraeder 28, 29, 30 ist in Fig. 3 
5 dargestellt. Eine alternative Zerlegung dazu ist in der Fig. 4 dargestellt. 

Eine Zerlegung einer Pyramide 30 in zwei Tetraeder 31, 32 mit den Eck-Knotenpunkten 33, 34, 35, 36, 37 ist in der 
Fig. 5 dargestellt. Eine alternative Zerlegung ist in der Fig. 6 gezeigt. 

In der Fig. 7 ist ein Tetraeder 40 dargestellt, der entweder Idealisierungs-Element oder ein Teilelement ist, in das bzw. 
in die Volumenelemente wie Hexaeder, Pentaeder, Pyramiden erfindungsgemaB zerlegt werden. Der Tetraeder hat die 
10 Eck-Knotenpunkte 41, 42, 43, 44 und die Flachen-Mittenknolen 46, 47, 48, 49, wobei der Flachenmitten-Knoten 49 in 
der Darstellung der Fig. 7 durch den Eck-Knotenpunkt 42 verdeckt und nicht gezeigt ist. 

Zur Durchfuhrung der erfindungsgemaBen Berechnungsverfahren anhand der Tetraeder- Volumenelemente (Fig. 7) 
werden die vier Eckpunkte der jeweiligen Tetraeder, d. h. die primaren Knotenpunkte, und sogenannte Flachenmitten- 
knoten, d. h. sekundare Knotenpunkte, als geometrische Bezugspunkte oder Strukturkoppelpunkte verwendet. Dabei 
15 werden erfindungsgemaB allgemcine Tetracder-Formcn, insbesonderc unregclmaBigc und auch schlankc, spitze Formcn 
verwendet. 

Allgemein basieren Finite-Elemente-Verfahren auf einem Energieansatz, nach dem die durch auBere Krafte auf die 
Struktur und den durch diese bewirkten Wee gebildete auBere Energie mit einer inneren Energie der Struktur gleichge- 
setzt wird. Dabei ist die innere Energie = E e^, wobei E das jeweilige elastizitatsmodul und e gleich der Dehnung als Ab- 

20 lcitung der Verschiebung an einem beliebigen Punkt des Volumenelements ist. Mil Hilfe dieses Ansatzes laBt sich auf 
den Verschicbungszustand beliebiger Punkte innerhalb der Struktur schlieBen. Um den Verschiebungszustand in einem 
beliebigen Strukturpunkt zu berechnen, werden allgemein Funktionen fur den Verschiebungszustand bei einem beliebi- 
gen Punkt innerhalb eines finiten Elements aufgestellt. Bei diesen Ansatzfunktionen werden Funktionen mit unbekann- 
ten Polynomkoeffizienten verwendet, die in eindeutiger Beziehung zu den Freiheitsgraden in den Knotenpunkten des je- 

25 weiligen Volumenelements stehen. Diese Relation zwischen den unbekannten Polynomkoeffizienten und den Freiheits- 
graden an den Knotenpunkten des jeweiligen Volumenelements, d. h. den Knotenfreiheitsgraden, wird bcstimmt durch 
die geometrische Form des verwendeten Volumenelements und die Art der Ansatzfunktion. 

Bei der erfindungsgemaBen Idealisierung der zu berechnenden Struktur ist die idealisierte Struktur gegebenenfalls 
nach einer Zerlegung der zunachst verwendeten Volumenelemente in Tetraeder gemaB Fig. 1 bis 6 vollstandig aus Tetra- 

30 edem aufgebaut. Dabei werden fur den Verschiebungszustand, d. h. fUr die Verse hiebungen im gesamten Bereich eines 
Tetraeders als Ansatzfunktion fur jede Koordinatenrichtung ein vollstandiges Polynom driller Ordnung, die sogenannte 
kubische Polynom- Form, verwendet. Die Ansatzfunktion unter Berucksichtigung der Diskretisierungsmatrix beinhaltet 
eine geometrische Beschreibung des betreffenden Tetraeders, insbesondere dessen Verschiebungszustandes. Sie wird ge- 
bildet mittels bekannter Methoden unter Verwendung der sogenannten Pascal' schen Funktionalschemen. In verkiirzter 

35 Schreibweise hat dieses Polynom fiir den Verschiebungsansatz foleende Form: 

R = LfJf- {q} = a t + a 2 x + a 3 y + a4Z + a 5 x 2 + a$xy + ayjr . . . + a 2 oz , wobei R eine Verschiebungskomponete in einer Ko- 
ordinaten-Richtung an einer beliebigen Stelle im Tttraeder-Element, LfJ eine Zeilenmatrix fiir die Poly norn-Ansatzfunk- 
tionen und {q} eine Spaltenmatrix fur die Polynom-Koeffizienten ist. 

Die voranstehende Polynom-Ansatzfunktion bezieht sich dabei lediglich auf eine Raumkoordinate; bei Betrachtung 

40 aller drei Koordinaten ergeben sich somit fur jedes Volumenelement sechzig Unbekannte. 

Die Koeffizienten ai bis a2o je Raumkoordinate bilden die zunachst unbekannten Polynom- Ansatzkoeffizienten. Bei 
einem Ansalz in Form des vollstandigen kubi schen Polynoms fiir das Tetraedervolumen, wie er erfindungsgemaB vorge- 
sehen ist, ergeben sich dabei immer fur jede Koordinate zwanzig unbekannte Polynom- Ansatzkoeffizienten. Diese drei- 
mal zwanzig Polynom-Ansatzkoeffizienten sind fur jedes Tetraederelement zu bestimmen. Dafur sind mathematisch 

45 dreimal zwanzig Zuordnungen oder Bestimmungsgleichungen erforderlich. 

ErfindungsgemaB werden fiir jedes Tetraederelement diese dreimal zwanzig Zuordnungen aus zwei Beitragen gebil- 
det: Zum einen wird jeder Knotenpunkt beschrieben durch einen Verschiebungsvektor sowie durch die drei Ableitungen 
dieses Raumvektors nach den drei verwendeten Koordinaten, zum andern werden fur jeden Tetraeder sekundare Knoten- 
punkte definiert. Diese sekundaren Knotenpunkte sind Rachenmitten- Knoten oder Schwerpunkte aller vier Seitenfla- 

50 chen jedes Tetraederelements. Fiir jeden dieser Mittenknoten werden lediglich die Freiheitsgrade der Verschiebungen 
eingebracht, nicht aber deren Ableitungen. Somit ergeben sich pro Tetraeder und pro Koordinatenrichtung 16 + 4 = 20 
Zuordnungen. Werden insgesamt alle drei Koordinatenrichtungen betrachtet, ergeben sich somit 60 Zuordnungen. Diese 
Zuordnungen ergeben sich nach ublichen mathemati schen Verfahren, die im Stand der Technik bekannt sind. In Form ei- 
ner mathemati schen Formulierung haben sie folgende Form: 

55 {Verschiebungs- Vektor fiir die primaren und sekundaren Knotenfreiheitsgrade) = [Formmatrix oder Diskretisierungs- 
matrix, die die inviduelle Tetraedergeometrie enthalt] { Vektor der Polynomkoeffizienten}. 

Da diese Zuordnung, die durch die Diskretisierungs- oder Form-Matrix beschrieben wird, eindeutig ist, kann diese 
Gruppe von Zuordnungen invertiert werden, d, h. nach den Vektor mil dem unbekannten Polynomkoeffizienten aufgelost 
werden, wobei die VerschiebungsgroBen als unbekannte RechengroBen eingefuhrt werden. Auf diese Weise sind diese 

60 Polynomkoeffizienten aufgrund von Matrizenoperationen bestimmbar. 

Die Bildung von Formmatrizen aus gegebenen geometrischen Elementen, sowie die Definition von Zuordnungen zwi- 
schen Freiheitsgraden der verwendeten geometrischen Elemente und den Unbekannten von Ansatzfunktionen beziiglich 
deren Knotenpunkten, sind allgemein nach dem Stand der Technik bekannt. ErfindungsgemaB werden jedoch bestimmte 
Finite-Elementtypen, namlich Tetraeder, zur Idealisierung von zu berechnenden Strukturen verwendet, fur deren innere 

65 Verschiebungszustande ganz speziell ein Ansatz mit vollstandigem kubischen Polynom und 20 Ansatzfuntionen gewahlt 
wird. Zur Definition der erforderlichen Zuordnungen zur Bestimmung der Polynomkoeffizienten werden erfindungsge- 
maB primare und sekundare Knoten definiert. 

Dabei sind die primaren Knoten die geometrischen Eckpunkte der Tetraeder. Diese primaren Knoten werden mittels 
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deren raumlichen Verschiebung sowie der Ableitungen dieser raumlichen Verschiebung nach den Koordinatenrichtun- 
gen erster Ordnung definiert. Die Freiheitsgrade in den primaren Knoten sind wie folgt definiert: 
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Die GroBen x, y, z stehen fiir die drei Raumkoordinaten und die GrtiBen u, v, w fUr die Komponenten des zugeordneten 
Verschiebungsvektors. 10 

Die sekundaren Knoten gehen dagegen nur mit ihren Verschiebungen in die Zuordnungen ein. Diese werden mit Ver- 
fahren nach dem Stand der Technik zur Vereinfachung des weiteren Ablaufs durch Kondensation und/oder Reduktion am 
Verbundelement, z. B. arn Hexaeder nach den Fig, 1 und 2, am Pentaeder nach den Fig. 3 und 4, an der Pyraminde nach 
den Fig. 5 und 6, eliminiert. Wenn der Tetraeder als Volumenelement verwendet wird, werden die Flachenmitten-Knoten 
durch Reduktion eliminiert. Dies geschieht dadurch, daB die Verschiebungen in den Flachenknoten durch die Verschie- 15 
bungen in den benachbarten Eckknoten beschrieben werden. Bei der Kondensation haiidelt es sich um eine exakte Eli- 
mination der Gleichungen fiir die Flachenmitten-Knoten durch entsprechende mathematische Operationen. Bei der Re- 
duktion werden die Verschiebungen in den Flachenmitten-Knoten durch die davon unabhangigen \ferschiebungen in den 
benachbarten Eckknoten ausgedriickt. Da die Kondensationsverfahren exakt sind, besteht grundsatzlich das Ziel, soweit 
wie moglich die Fl ache n-Mittcn knoten durch Kondensation zu eliminiercn. Diese Verfahrensschritte sind Stand der 20 
Technik. 

In den nachsten Verfahrenschritten werden nach iiblicherweise verwendeten Energieprinzipien aus den ermittelten 
Gleichungen fur die Tetraeder-Eckknoten, die den Verschiebungszustand fiir den jeweiligen gegebenen falls aus dem ent- 
sprechenden Verbundelement gebildeten Tetraeder beschreiben, Volumenintegrale uber jeden dieser Tetraeder gebildet. 
Daraus werden die sogenannten Element-Streifigkeitsmatrizen fur jeden dieser Tetraeder definiert. Die Energieanteile je- 25 
des Tetraeders werden gcgebenenfalls zu den Energieanteilen eincs Verbundelements, also z. B. eines Hexaeders, eines 
Pentaeders oder einer Pyramide, zusammengesetzt. Aus den Verbundelement-Streifigkeitsmatrizen werden wiederum 
die System-Streifigkeitsmatrizen zusammengesetzt, wobei die gleichen Energiemethoden verwendet werden. Die Sy- 
stem-Streifigkeitsmatrizen bilden bekanntermaBen die Grundlage fiir die Losungsgleichungen des jeweils gestellten 
Strukturproblems. Beispielsweise k6nnen bei Vorgabe auBerer Krafte mittels der System-Streifigkeitsmatrix die Ver- 30 
schiebungen in den Knotenpunkten der idealisierten Struktur berechnet werden. 

Zur Verbesserung der Ergebnisgenauigkeit kann zusatzlich folgender Verfahrensschritt vorgesehen werden: 
An Stellen der zu berechnenden Struktur an denen diese mittels Verbundelementen idealisiert worden sind, die in Tetra- 
eder als Teilelemente zerlegt werden, erfolgt diese Zerlegung erfindungsgemaB je Verbundelement aufgrund beider Zer- 
legungsvarianten, Beispielsweise erfolgt die Zerlegung eines zur Idealisierung verwendeten Hexaeders nach der ersten 35 
(Fig. 1) und in einer weiteren Idealsierung nach der zweiten Variante (Fig. 2). Ebenso erfolgt die Zerlegung jedes Panta- 
eders sowohl nach Fig. 2 als auch nach Fig. 3. Die Zerlegung der Pyramide erfolgt sowohl nach Fig. 5 als auch nach Fig. 
6. Dadurch wird eine durch den Benutzer idealisierte Struktur durch das erflndungsgemaBe Verfahren auf zwei sich er- 
ganzende Weisen in Tetraeder zerlegt. Auf diesen beiden aufgrund verschiedener Zerlegungen der verwendeten \fer- 
bundelemente hervorgegangenen Teilelement-Strukturen werden die beschriebenen Berechnungs verfahren jeweils ei- 40 
gens angewendet, so daB jede idealisierte Struktur zweimal berechnet wird. 

Auf diese Weise wird jeder Knotenpunkt der idealisierten Struktur systematisch gleich gewichtct, da jede Vierecksfla- 
che eines Verbundelements auf beide moglichen Arten in Dreiecksflachen zerlegt wird. Bei einer Zerlegung eines \fer- 
bundelements, z. B. eines Hexaeders, sind einzelne Eck-Knoten einer viereckigen Seitenftache Teilnehmer an mehreren 
Tetraeder-Elementen. Z. B. ist der Eck-Knoten 2 des Hexaeders der Fig. 1 an den Tetraedem 11, 13, 15 beteiligt, wahrend 45 
der Eck-Knoten 3 nur dem Tetraeder 11 zugeordnet ist. Werden die beschriebene Berechnung nur aufgrund einer be- 
stimmten Zerlegung eines Verbundelements durchgefiihrt, ergibt sich numerisch eine unterschiedliche Bewertung von zu 
mehreren Tetraedern gehorende Eck-Knoten im Vergleich zu anderen. Diese unterschiedliche Bewertung fuhrt zu Unge- 
nauigkeiten, die durch die zweifache Idealisierung beseitigt werden konnen. 

Berechnungstechnisch werden bei der Bildung von Teilelement-Strukturen nach beiden Varianten die jeweils ermittel- 50 
ten Teilelemente (Tetraeder)-Steifigkeiten addiert. Die Steifigkeiten der Tetraeder, die zur Idealisierung der Struktur ur- 
spriinglich verwendet worden sind, werden verdoppell. Daraus ergibt sich in etwa der doppelte Wert fiir die jeweilige 
Teilelement-Steifigkcit bzw. deren inneren Energie sowie fiir die Gesamt-Steifigkeit bzw. die innere Energie der gesam- 
ten Struktur. Dieser doppelte Wert wird im weiteren Verfahren dadurch eliminiert, daB beispielsweise dieser Wert wie- 
derum halbiert wird oder die auBeren Lasten verdoppelt werden. Auch andere mathematische Verfahren sind moglich. 55 

Das erflndungsgemaBe Verfahren wird im folgenden zusammenfassend anhand der Fig. 8 beschrieben: Bei dem erfin- 
dungsgemaBen Verfahren zur Berechnung von Strukturen mit Hilfe der Finiten-Elemente-Methode wird vorzugsweise 
zunachst eine Finiten-Elemente-Methode (FEM)-Vorverarbcitung 101 durchgefiihrt. Dabei wird die zu berechnende 
Struktur fur nachfolgende Schritte vorbereitet, wobei z. B. Konturlinien oder AuBenkonturen geglattet oder unwichtige 
Details der zu berechnenden Struktur entfernt werden. Das durch die FEM- Vorverarbeitung 101 aufbereitete Modell 102 60 
wird dann in einem Schritt der Idealisierung 103 durch ein Maschennetz beschrieben. Das Maschennetz definiert sich 
durch die Idealisierung der zu berechnenden Struktur durch Volumenelemente, die aus Hexaedern, Pentaedem, Pyrami- 
den oder Tetraedern gestaltet sind. Die GroBe und Form der Vqlumenelemente ergibt sich aus der Form der zu berech- 
nenden Struktur bzw. des aufbereiteten Modells 102, wobei das aus der Idealisierung 103 entstehende idealisierte Modell 
104 mit seinen Eck- Knotenpunkten von der Kontur des aufbereiteten Modells 102 bestimmt ist. Das idealisierte Modell 65 
104 stellt eine mbglichst gute Annaherung an die raumliche Gestalt des aufbereiteten Modells 102 dar. Es konnen bei der 
Idealisierung jedoch auch andere Aspekte nach bestehendem Fachwissen eine Rolle spielen, z. B. angenommene Kraft- 
verlaufe oder Gesichtspunkte der Lagerung der Struktur. Das idealisierte Modell 104 ist - soweit es das beschriebene Be- 
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rechnungsverfahren betriffit - ausschlieBlich aus Volumenelemente gebildet. Dabei haben die Volumenelemente, die 
selbst wieder aus Tetraedern zusammengestzt werden konnen, in den nachfolgenden Schritten die Funktion von Ver- 
bundelementen. In einem nachfolgenden Schritt erfolgt eine Zedegung 105 der Volumenelemente, soweit sie in dem 
idealisierten Modell 104 nicht als Tetraeder ausgestaltet sind, bzw. der Verbundelemente in Tetraeder nach den voranste- 

5 hcndcn anhand der Fig. 1 bis 6 gcmachten Ausfuhrungen. Das Ergcbnis der Zcrlcgung 105 ist, das idcalisierte Modell 
104 in eine Teil-Element-Struktur 106 ausschlieBlich aus Tetraedern umgeformt ist. Gegebenenfalls werden Zusatzinfor- 
mationen 107 iiber die verwendeten Volumenelementen, also iiber die verwendeten Hexaeder, Pentaeder und Pyramiden, 
beriicksichtigt und der Zerlegung 105 zugefiihrt. In einem nachfolgenden Schritt 108 erfolgt die Beschreibung des Ver- 
schiebungszustands fUr jedes Tetraeder-Element durch ein vollstandiges Polynom 109 dritter Ordnung je Koordinanten- 

10 richtung. In einem weiteren Schritt 110 erfolgt die Bestirnmung der Polynom-Koeffizienten durch Zuordnungen zwi- 
schen den Knotenfreiheitsgraden, Faktoren zur Beschreibung der individuellen Tetraeder-Geometrie und dem jeweiligen 
Polynom-Koeffizienten. Als Knotenfreiheitsgrad wird fiir jeden Knoten-Eckpunkt der Verschiebungsvektor und die drei 
Ableitungen dieses - Vektors sowie der Verschiebungsvektor fUr die Flachenmitten-Knoten jedes Tetraeder-Elements ver- 
wendet. In einem Schritt 111 werden die Teilelement-Tetraeder, soweit sie im Schritt 105 aus einer Zerlegung hervorge- 

15 gangen sind, wicdcr in die Volumenelemente, wie sie im Schritt 105 erzeugt worden sind, bzw. in Verbundelemente zu- 
sammengesetzt, so daB sich ein Verbundaufbau 112 daraus ergibt. Gleichfalls wird im Schritt 111 eine Elimination der 
Flachenmitten-Knotenpunkte vorgenommen. Im nachfolgenden Schritt 113 erfolgt der Aufbau der Systemmatrix aus den 
Verbundelementen 112, d. h. die Verschiebungen und Ableitungen werden iiber alle Volumenelemente aufsummiert und 
dadurch die Systemmatrix gebildet. AbschlieBend erfolgt im Schritt 115 die Losung des gestellten mechanischen Pro- 

20 blems mittels der Systemmatrix 114. 

In dem bestehenden Verfahren ist eine Tetraeder-Schleife 121 und cine Volumenelemcnt-Schleife 122 eingebaut. Bei 
der Tetraeder-Schleife 121 geht es um die Berechnung der Tetraeder, aus denen die Volumenelemente zusammengesetzt 
sind. Die Tetraeder-Schleife 121 ist also nur erforderlich fur die Volumenelemente, die im Schritt 105 in Tetraeder zerlegt 
worden sind. Fur jeden aus einem Volumenelement in der Zerlegung 105 entstehenden Tetraeder ist eine Schleife 121 

25 durchzufuhren, Auf diese Weise liegen als Ergebnis des Schritts 111 die Steifigkeiten fur die Volumenelemente, wie sie 
im idealisierten Modell 104 vorgesehen warcn, vor. Um im Schritt 113 die Systemmatrix zu bilden, ist die Volumenele- 
ment-Schleife 122 erforderlich. Es ist fur jedes Volumenelement, wie es im idealisierten Modell 104 vorgesehen ist, ein 
Anteil an Steifigkeit oder Energie fur die Systemmatrix vorzusehen und je Volumenelemt aufzuaddieren. Nachdem die 
Schleife 122 fur jedes Volumenelement durchlaufen ist, sind alle Anteile in der Systemmatrix 114 eingebracht. Der ab- 

30 schlieBende Schritt der Systemlosung und der Nachbearbeitung 115 erfolgt nach Verfahren nach dem Stand der Technik, 
Zur Systemlosung gehort die Berechnung der unbekannten Verschiebungen und zur Nachbearbeitung gehort die Aus- 
wertung der Verschiebungs- und Spannungsergebnisse. 

Das erfindungsgemaBe Verfahren kann mit anderen Verfahren nach dem Stand der Technik kombiniert sein. So kann 
ein Teil einer zu berechnenden Struktur mit dem erfindungsgemaBen Verfahren und ein anderer Teil der Struktur mit ei- 

35 nem Verfahren nach dem Stand der Technik berechnet werden, bei dem auch andere Elemente als Volumenelemente zur 
Idealisierung der Struktur verwendet werden. 

Ein Vorteil der Erfindung ist, daB die kleinsten Volumenelemente, iiber dercn Volumcn zur Bildung dercn Steifigkeits- 
matrizen integriert wird, erfindungsgemaB ausschlieBlich Tetraeder sind, so daB sich eine hohe numerische Stabilitat des 
erfindungsgemaBen Verfahrens ergibt, da Tetraeder sogenannte Simplexformen sind, iiber die mathematisch exakt, bei- 

40 spielsweise mittels Bessel-Funktionen, integriert werden kann. Nur bei Verwendung von Tetraedern als Idealisierungs- 
oder auch als Teil-Elemente, konnen die Polynom- Ansatzfunktionen dritter Ordnung numerisch stabil verarbeitet wer- 
den. Zur Eliminierung der Flachenmittenknoten der Tetraeder-elemente werden wahlweise Kondensations- oder Reduk- 
tionsverfahren nach dem Stand der Technik angewendet. Vorzugs weise werden Kondensationsverfahren angewendet, 
die exakt sind, so daB sie die numerische Stabilitat des Programms verbessern. 

45 Die Verwendung von Tetraedern als kleinste Volumenelemente, iiber die zu den Elementsteifigkeitsmatrizen integriert 
werden, liegt darin, daB auch sehr schlanke Tetraeder, d. h. Tetraeder mit einem Hohen-Seitenverhaltnis von iiber 100, 
noch numerisch stabil sind. Der Grund liegt darin, daB die Tetraeder numerisch exakt integriert werden konnen. 

Durch die hohe Rechengenauigkeit, die das erfindungsgemaBe Verfahren auf der Basis der vollstandigen kubischen 
Polynom-Ansatzfunktionen unter Einfuhrung der vollen Ableitungsfreiheitsgrade liefert, wird die erfolgreiche Analyse 

50 aller bekannten statischen und dynamischen Festigkeits- und Stabilitatsprobleme ermoglicht. Unter letzterem Begriff 
sind auch die mathematisch sensiblen Formulierungen zur geometrisch nichtlinearen Problematik der groBen \ferfor- 
mungen und der Strukturstabilitat zu verstehen. 

Die hohe Ergebnisqualitat des erfindungsgemaBen Verfahrens aufgrund der kubischen Polynom-Ansatzfunktionen er- 
offnet die Moglichkeit, alle auf der Basis von 3D-Modellen nach heutigen Verfahren erstellte Konstruktionen direkt, d. h. 

55 mit einem sehr geringen Umformungs- und Formalisierungsaufwand, als 3D- Volumen-Modelle in der FE-Berechnungs- 
technik zu simulieren. Die bislang ubliche Modellierung der Strukturen als Stabwerke, spezielle Schalen, Faltwerke und 
3D-Teilanalysen kann entf alien und wird durch eine einheitliche Volumenmodellierung ersetzt. Durch diese Vorgehens- 
weise werden anwenderabhangige und modellierungsbedingteUngenauigkeiten weitgehend reduziert. Der hohe Ansatz- 
grad der Polynom-Ansastzfunktionen ermoglicht es, relativ grobe Maschennetze zur Idealisierung der zu berechnenden 

60 Struktur einzufuhren. So geniigt iiber die Dicke allgemeiner Schalenstrukturen eine Elementreihe aus Verbundelemen- 
ten, z. B. Hexaedem, zur Idealisierung der Struktur. 

Die Spannungen in den Knotenpunkten der idealisierten Struktur, werden durch das erfindungsgemaBe Verfahren ein- 
deutig berechnet und brauchen nicht wie bei konventionellen Verfahren elementabhangig ermittelt zu werden. Der Grund 
liegt darin, daB bei dem erfindungsgemaBen Verfahren die Ableitungen der Verschiebungen in den Systemknoten der 

65 idealisierten Struktur eindeutig definiert sind, so daB daraus auch die Dehnungen und iiber die Werkstoffgesetze auch die 
Spannungen in diesen Punkten eindeutig definiert sind. Da diese Losungen in bezug auf die System-Knotenpunkte vor- 
liegen, sind sie unabhangig von dem zu dem jeweiligen System-Knotenpunkt gehorenden Finite-Element. Dadurch er- 
iibrigt sich die Mittelwertbildung bzw. Extrapolation an diesen System-Knotenpunkten, wie sie im Stand der Technik 
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vorgenommen werden muB, da die Dehnungen, Verschiebungen und Spannungen an diesen Punkten elementabhangig 
ausgerechnet werden. Dies ist dadurch moglich, da die Ableitungen in den Freiheitsgraden verwendet werden und in die 
Polynom-Ansatzfunktionen eingegangen sind, rechentechnisch direkt erfafit werden konnen. Demgegenuber sind beim 
Stand der Technik die Ableitungen der Verschiebungen naherungsweise ftlr die Elemente auf der Systemebene in einem 
Postprozcssingverfahren zu bcstimmen. 5 

Die Brauchbarkeit der erfindungsgemaB ermittelten Spannungen an dem System-Knotenpunkten wird auch dann ge- 
wahrleistet, wenn bei der Idealisierung der Struktur ein grobes Maschennetz verwendet wurde, da die Spannungen nicht 
aus einer Mittelwertbildung oder Extrapolation ermittelt werden. - Daraus ergeben sich signifikante Vorteile bei Analysen 
zu einem RiBfortschritt oder zu Kerb spannungen in der zu berechnenden Struktur, weil Spannungsspitzen erflndungsge- 
maB zuverlassig identifiziert werden. 10 

Dariiber hinaus besteht durch das erfindungsgemaBe Verfahren die Moglichkeit, die Automatisierung fur die Erstel- 
lung des Finite-Element-Modells vollstandig zu perfektionieren und weitgehend unabhangig von den individuellen Ge- 
pflogenheiten des Anwenders zu machen. Neben der Steigerung der Ergebnisgenauigkeit werden dadurch potentionelle 
Fehlerquellen durch unsachgemaBe Struktur-Idealisierungen ausgeschiossen. 

Bei Verfahren nach dem Stand der Technik vernachlassigte gcomclrische Strukturdetails wic Dickcnspriinge, Exzcnt- 15 
zitaten, Radiusiibergange konnen in der Finite-Element-Methode wegen der Idealisierung der Struktur ausschlieBlich mit 
Volumenelementen quantitativ beriicksichtigt werden. Auch Materialinhomogenitaten, wie diese bei Faserverbundwerk- 
stoffen und Plasuzierungen auftreten, konnen simuliert werden. 

Patentanspruche 20 

1. Verfahren zur Berechnung von Strukturen mit Hilfe der Finiten-Elemente-Methode mit den Schritten: 

- Idealisierung der zu berechnenden Struktur durch Voiumenelemente, die die Form eines Hexaeders, Penta- 
eders, einer Pyramide oder eines Tetraeders besitzen, und Festlegung der Eck-Knotenpunkte in der idealisier- 

ten Struktur, 25 

- Zerlegung der Voiumenelemente mit Hexaeder-, Pentaeder- und Pyramiden-Form in Tetraeder, so daB die zu 
idealisierende Struktur nur Teilelemente in Tetraeder-Form aufweist, 

- Beschreibung des Verschiebungszustands fur jedes Tetraeder-Element durch ein vollstandiges Polynom 
dritter Ordnung je Koordinatenrichtung, 

- Bestimmung der Polynom-Koeffizienten durch Zuordnungen zwischen den Knotenfreiheitsgraden, Fakto- 30 
ren zur Beschreibung der individuellen Tetraedergeometrie und den jeweiligen Polynom- Koeffizienten, wobei 

als Knotenfreiheitsgrad fur jeden Knoten-Eckpunkt der Verschiebungsvektor und die drei Ableitungen dieses 
Vektors, und der Verschiebungsvektor fur die Flachenmitten-Knoten jedes Tetraederelements verwendet wird, 

- Elimination der Flachenmitten-Knotenpunkte, 

- Addition der Verschiebungen und Ableitungen jedes letraeder-Teilelements und Bildung der Systemmatrix, 35 

- Losung des gestellten mechanischen Problems mittels der Systemmatrix. 

2. Verfahren zur Berechnung von Strukturen mit Hilfe der Finiten-Elemente-Methode nach dem Anspruch 1 , da- 
durch gekennzeichnet, daB die Ansatzfunktionen dritter Ordnung vollstandig sind und mittels des Pascal'schen 
Funktionalsschema erstellt werden. 

3. Verfahren zur Berechnung von Strukturen mit Hilfe der Finiten-Elemente-Methode nach den Anspruchen 1 oder 40 
2, dadurch gekennzeichnet, daB die Elimination der Flachenmitten-Knoten mittels Reduktions verfahren erfolgt. 

4. Verfahren zur Berechnung von Strukturen mit Hilfe der Finiten-Elemente-Methode nach einem der voranste hen- 
den Anspriiche, dadurch gekennzeichnet, daB die Elimination der Flachenmitten-Knoten mittels Kondensations ver- 
fahren erfolgt. 

5. Verfahren zur Berechnung von Strukturen mit Hilfe der Finiten-Elemente-Methode nach einem der voranstehen- 45 
den Anspriiche, dadurch gekennzeichnet, daB die Bildung der Teilelemente aus einem Volumenelement in zwei sich 
erganzenden geometrischen Varianten erfolgt und die nachfolgenden Schritte aufgrund einer Uberlagerung der aus 
beiden Varianten resultierenden Ergebnisse erfolgt. 
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The invention relates to a method of calculating structures using the finite element method. The inventive 
method comprises the following steps: idealizing the structure to be calculated by volume elements and 
fixing the corner nodes in the idealized structure, breaking down the volume elements that have a 
hexahedron, pentahedron or pyramidal shape into tetrahedrons so that the structure to be idealized is a 
partial element structure made up of tetrahedrons, describing the displacement state for every tetrahedron 
element by a complete third order polynomial per coordinate axis, determining the polynomial coefficients 
by establishing correlations between the node degrees of freedom, the factors for describing the individual 
tetrahedron geometries and the respective polynomial coefficients, with the displacement vector and the 
three derivatives of said vector being used as the node degree of freedom of every node corner and the 
displacement vector for the centroid nodes of every tetrahedron element being used, eliminating said 
centroid nodes, integrating and solving the respective mechanical problem by means of the system 
matrix. 
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